1  .. 

I  AD-A106  37 

n  UNCLASSIFZ 

i 

ED 

AIR  FORCE  INST  OF  TECH  WRI6HT-PATTERS0N 
COMPARISON  OF  THE  DYNAMIC  PROCESSES  DURI 
JUN  81  J  H  DAVENPORT 

AFIT-CI-81-29T 

AFB  OH 

NS  MAJOR  AND 

F/G  4/1 

"INOR  STRA— ETC(U) 

NL 

m 

F 

1 

■■ 

i 

— 

-J 

l 

1 

1 

■■■ 

■■■ 

■■1 

■1 

■  _ 

ME  nil  copy  ad  a  10*;  :m 


81-?9T 


*  !  «  S4  Milk  ASM  SSHUr 


’*•  IA'»  ♦.m  -••■»'«  '  «  '  n.n'i  '  •  »*«  ,  ,mii  ■  tx.  t  »  .  *  irvri'.' 

*.  .  .«*  1  » ta*t  t  .  -I  t  .  ■  •»,,.<>  .  ♦  ( •**  A '  •  *  *  f  !*»  *  >  t  *t«  *  tw*i  ,N,  A*,  :  t  •,,,<)>* 

A*  »»t  .  *  ’  *  !  MJ  ’  .  .  •  H,  .1  •  '  «  •  <>H  \  »  .«>«>>>  Jfl.j  •  t  . 


**  S» 

•*  '  i*  t  <  A  1  I  #•  •  V  i  -«  A'  *  ,  *•  t  '  * 

Hiin  i't  the  i>ynj»n  Prtxrssrs  During  Haior  and  Hi*>or  Statniptwrn 

Umin^i  i"  A  HumpiuaI  Nod*  1 

a,  «,»  -1a*»s  fa«*npo> t 

>t'|M  ■  4"  t  S  ■  *1' 


tw  t  *  #  *  ' 


9  O*  •  ‘  $r  1  ?  i  »  .,g»  *  **'  ‘  tr 


mw  * 


*«••♦*  *  *  A  .  «  .  i  .« v  * 

*  nfi  *  «  *  r » f**1  »  *  .  r ■<*  h ,  v  1  ’  *  4 

**•*»..►  «*v  ' -4  *•#%•  *  1  *  *  a-*  ?■#*« 

*  i  *  *4»A,  >*r*  V*  »1«  A' 


■  •  A*- 


4  .>•»«» 

’  fVw  tv  »'«  Ihf 


A*  w*  v’* 
4*t  A  1  v  C»*»  *rr  ns  *  *  • 
r  ■  ,<*r  th#  • 4  r 


’  N’  * 


r  *  T  .  <*  it  f  ►  •  ,  i  ,  > 

,  s  J  '  II  A  f  ♦  ’  *  *  V 

►  ««  f.fif  i  1,*.#  hf‘..s# 


*  »  *  •  * 
w  •  « 

♦  i  ‘  /•  •  i 


•aN 


.  .  s»H  *  ' 

,  *N  ftN 


*»  s 

A  ,  \\ 


>>*  * «  .»*■  *  **  M>v*  i4»f  v 

,  *  . f  ♦  ;>ot  #>»  *  •  a  >  *  s  a 

t  '  *  •  r  »  ,',J'  '  *  •  *  nwir*  * 


■KAT't 

\  V  •  >, 


V  •  *  '  V 


f!M[S' 


Q'  Vl  I 


\ 


COMPARISON  Of  THE  DYNAMIC  PROCESSES  DURING 
MAJOR  ANO  MINOR  STRATOSPHERIC  WARMINGS 
IN  A  NUMERICAL  MOOEl 


t>y 

Jamas  Howard  Davanport 


An  abstract  of  a  thasls  subwit tad  to  tha  faculty  of 
Tha  Uni  varsity  of  Utah 

In  partial  fulflllmant  of  tha  raqulramants  for  tha  dagraa  of 


Mastar  of  Sclanca 


S.  k.  «cao 


Chairman,  Suparvlsory  Commlttae 
Profassor  of  Nataorology 


Dapartmant  of  Mataorology 
Tha  Uni  varsity  of  Utah 
Juna  IRfll 

8  i  o 


J 


*  *  \  o 


ABSTRACT 


In  this  study  we  have  analyzed  the  processes  Involved  In  a  major 
and  minor  warming  based  on  data  generated  from  a  31-layer  primitive 
equation  spectral  model  consisting  of  a  5-layer  troposphere,  12-layer 
stratosphere  and  14-layer  mesosphere.  The  forcing  of  the  model  con¬ 
sisted  of  orographic  features  designed  to  simulate  the  major  Northern 
Hemisphere  continents.  This  gave  a  combination  of  wave  number  1  and 
wave  number  2  forcing.  The  analysis  of  the  simulated  major  and  minor 
warmings  Indicated  several  similarities  and  differences  between  the 
cases. 

The  Initial  difference  between  the  cases  Is  that  case  1  (major 
warming)  had  a  weak  polar  vortex,  while  case  2  (minor  warming)  had  a 
strong  polar  vortex.  This  resulted  In  case  2  having  much  larger 
turbulent  flux  terms,  but  not  sufficient  to  reverse  the  mean  polar 
winds.  Case  1  produces  a  major  warming  In  the  polar  stratosphere. 

Similarities  of  the  cases  Included:  wave  number  1  waves 
dominant  In  the  mesosphere  and  wave  number  2  waves  dominant  In  the 
stratosphere;  development  of  a  meridional  circulation  of  an  approxi¬ 
mate  8  to  12  day  cycle  In  the  stratosphere  and  mesosphere,  and  a  good 
correlation  between  the  peaks  of  this  meridional  circulation  and 
change  In  the  polar  night  jet.  The  peaks  of  the  meridional  circula¬ 
tion  coincided  with  magnitude  peaks  of  convergence  and/or  divergence 
of  momentum  and  vertical  flux  of  geopotentlal . 

An  analysis  of  the  latitude-height  sections  of  meridional  flux 
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ABSTRACT 


In  this  study  we  have  analyzed  the  processes  Involved  In  a  major 
and  minor  warming  based  on  data  generated  from  a  31-layer  primitive 
equation  spectral  model  consisting  of  a  5-layer  troposphere,  12-layer 
stratosphere  and  14-layer  mesosphere.  The  forcing  of  the  model  con¬ 
sisted  of  orographic  features  designed  to  simulate  the  major  Northern 
Hemisphere  continents.  This  gave  a  combination  of  wave  number  1  and 
wave  number  2  forcing.  Ihe  analysis  of  the  simulated  major  and  minor 
warmings  indicated  several  similarities  and  differences  between  the 
cases. 

The  Initial  deference  between  the  cases  Is  that  case  '.  major 
warming)  had  a  weak  polar  vortex,  while  case  2  (minor  wanning'  had  a 
strong  polar  vortex.  This  resulted  in  case  2  having  much  larger 
turbulent  flux  terms,  but  not  sufficient  to  reverse  the  mean  pc'ar 
winds.  Case  1  produces  a  major  warming  In  the  polar  stratosphere. 

Similarities  of  the  cases  Included:  wave  number  1  waves 
dominant  In  the  mesosphere  and  wave  number  2  waves  dominant  in  the 
stratosphere;  development  of  a  meridional  circulation  of  an  approxi¬ 
mate  8  to  12  day  cycle  In  the  stratosphere  and  mesosphere,  and  a  good 
correlation  between  the  peaks  of  this  meridional  circulation  and 
change  In  the  polar  night  jet.  The  peaks  of  the  meridional  circula¬ 
tion  coincided  with  magnitude  peaks  of  convergence  and/or  divergence 
of  momentum  and  vertical  flux  of  geopotential. 


An  analysis  of  the  latitude-height  sections  of  meridional  flu* 


of  sensible  heat  and  vertical  of  geopotentlal  Indicated  that  the 
relationship  Between  vertical  flux  of  geopotentlal  and  meridional 
flux  of  sensible  heat,  derived  from  the  linearized  equations,  does 
not  hold  for  the  nonstationary  developing  stage  of  the  stratospheric 
warming.  This  Inconslstancy  casts  questions  on  the  applicability 
of  the  concept  of  the  critical  level  to  the  state  of  nonlinear  Inter¬ 
action. 
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CHAPTER  1 


INTRODUCTION 

Records  of  winter  stratospheric  warmings  Indicate  that  about 
half  of  the  warmings  are  classified  as  major  warmings  defined  as  the 
occurrence  of  changing  tonal  westerly  to  easterly,  the  other  half 
are  classified  as  minor  warmings  (Mclnturff,  1978).  Analyses  of 
major  and  minor  stratospheric  warmings  show  that  during  major  warm¬ 
ings  both  the  temperature  and  tonal  wind  reversals  occur  first  in  the 
upper  stratosphere  then  extend  downward  to  the  lower  stratosphere  and 
troposphere,  whereas  during  minor  warmings  no  reversals  of  tempera¬ 
ture  and  tonal  wind  occurs  In  the  stratosphere,  and  the  tropospheric 
westerly  Jet  remains  almost  unchanged.  Indicating  little  feedback 
from  the  stratosphere  to  the  troposphere  (Koermer  and  Kao,  1980). 

In  view  of  the  distinctly  different  stratosphere-troposphere  Interac¬ 
tions  during  the  stratospheric  major  and  minor  warmings,  a  numerical 
simulation  of  the  stratospheric  major  and  minor  warmings  with  the  use 
of  31 -level  primitive  equations  spectral  model  has  recently  been 
performed  (Koermer,  1980).  It  Is  found  that  results  obtained  by 
Integrating  the  numerical  model  agree  well  with  those  of  the  observed, 
and  that  nonlinear  Interaction  plays  an  Important  role  in  the  develop¬ 
ment  of  stratospheric  major  and  minor  warmings.  The  objectives  of 
this  study  are  to  analyze  the  numerical  model  by  comparing  two  cases 
from  it.  The  first  case  Involved  a  major  stratospheric  warming.  The 
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second  case  Involved  a  minor  stratospheric  Manning.  By  examining 
these  numerical  simulations,  we  hope  to  gain  a  better  understanding 
of  the  dynamic  processes  Involved  In  the  evolution  of  the  strato¬ 
spheric  major  and  minor  warmings.  Specifically,  we  will  analyze  (1) 
the  effect  of  vertical  flux  of  geopotentlal  generated  from  orographic 
forcing  on  the  development  of  wave  motion  and  meridional  fluxes  of 
zonal  momentum  and  sensible  heat  In  the  stratosphere  and  mesosphere; 
(2)  the  development  of  zonal  mean  meridional  circulation,  as  the 
consequence  of  the  convergence  of  meridional  fluxes  of  zonal  momentum 
and  sensible  heat;  (3)  the  evolution  of  the  zonal  mean  maxima  of 
zonal  velocity,  temperature  and  variances  and  covariances  of  velocity 
components  and  temperatures;  and  (4)  the  effects  of  waves  of  various 
wave  numbers  on  the  development  of  major  and  minor  stratospheric 
warmings. 

The  focus  of  our  comparison  Is  primarily  on  the  zonal  mean  of  u, 
v,  w,  T,  «,  deviations  from  that  zonal  mean,  and  cross  correlations 
of  these  quantities  form  day  10  to  day  40  of  each  case.  Throughout 
this  paper  the  zonal  mean  eddy  flux  quantities  will  be  referred  to 
simply  as  zonal  mean  fluxes. 
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CHAPTER  2 


THE  PRIMITIVE  EQUATION  SPECTRAL  MOOEL 

The  model  used  In  this  study  basically  consists  of  two  parts. 
The  upper  part  of  the  model  Includes  the  stratosphere  and  mesosphere 
In  which  a  log-pressure  coordinate  system  Is  used.  The  lower  part 
of  the  model  represents  the  troposphere  In  which  a  modified  sigma 
coordinate  system  Is  used  In  order  to  handle  orography  of  the  earth. 

Figure  1  depicts  the  vertical  grid  structure  of  the  coupled 
models.  In  the  stratosphere  and  mesosphere,  az  *  3km  with  26  prog¬ 
nostic  levels,  which  will  generally  be  represented  with  the  Index  K . 
The  Index  r  Is  used  to  Indicate  the  five  prognostic  levels  of  the 
troposphere  where  Aa  •  .2.  Vertical  velocity  (w  or  a )  and  geopo- 
tentlal  (*)  are  diagnostic  variables  and  values  are  computed  for 
diagnostic  levels  centered  between  the  prognostic  levels.  Vortlclty 
(c),  divergence  (0)  and  temperature  (T)  are  the  prognostic  variables 
common  to  both  models.  Additionally,  log  pressure  (q)  Is  a  prog¬ 
nostic  variable  In  the  troposphere.  For  boundary  conditions,  we 
assume  vertical  velocity  Is  zero  at  the  top  of  the  upper  model  and 
the  bottom  of  the  lower  model. 

The  large  scale  dynamics  of  an  atmosphere  are  basically 
governed  by  the  equations  of  horizontal  motions,  thermocjynamlcs, 
mass  continuity,  and  hydrostatic  equilibrium.  For  the  tropospheric 
part  of  the  model,  the  governing  equations  In  the  sigma  coordinate 
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system  may  be  written  as  follows: 


dV--fkxv-7a-  gRT  7P„  ♦  F  , 
at  oP  +p 

*  R 


dT  .  RT  d  *  Q 
71  c  (o'P  *p  )  c" 

p  *  R  p 


(1) 


(2) 


3P 

*  ♦  7  •  IP  y)  +  P  35  »  0  i  (3) 

TF  *  *  Jo 


where 


li 

3o 


RTP 

★ 

aP  +“p“  ’ 
*  R 


o  -  (p-PR)/{p$-pR)  • 


(4) 


(5) 


P  is  a  reference  pressure  Indicating  the  upper  boundary  of  the  tropo- 

K 

spheric  model  and  lower  boundary  of  the  upper  atmospheric  model; 
p  Is  the  surface  pressure;  P*  *  p.  -  PD  ;  v  Is  the  horizontal 

w  j  K 

velocity;  f  Is  the  Coriolis  parameter;  k  Is  the  vertical  unit 

♦ 

vector;  T  Is  the  temperature;  F  Is  the  frictional  force;  Q  Is  the 

diabetic  heating;  c  Is  the  specific  heat  of  dry  air  at  constant 
P 

pressure;  p  *  dp/dt  Is  the  vertical  velocity  In  pressure  coordinates; 
o  ■  do/dt  Is  the  vertical  velocity  In  o-coordl nates;  t  represents 
time;  and  7  is  the  horizontal  gradient  operator. 

For  the  upper  atmosphere,  the  governing  equations  In  the  log- 
pressure  coordinate  system  may  be  written  as  follows: 


L 


■  mi/tti  viai— 


^  -  i  -  A"  J  *. 


where 


,  ♦  +  ♦ 
jJv*-fkxv-V$  +  F, 
3t 


dT  .  Q  .  N  wH 

Si  T  “R  * 

P 


7  .  y  +  3w  _  w  »o* 

32  TT 


d*  RT 
3«  .  RT  0  ,  o 
a  I  IT  *  ST  “FT 


d  s  3 _ +  v  •  v  +  w  ® _  » 

<IF  at  3z 


2  »  -  H*n  (p/p  )  , 

O 


w  *  d*  » 

cfF 


H  •  RT  /g  , 
$ 


2  dT  <7 
N  =  R  (  o  +  0) 

W  W  IT 


K  5  R/C 


3y  the  Hemholt2  theorem  v  can  be  expressed  as  the  sum  of  non 


divergent  and  Irrotatlonal  parts  as  follows 


V  «  k  X  V\f»  ♦  vx  , 


where  Is  the  stream  function  and  x  is  the  velocity  potential. 
Hence,  It  follows  that 
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♦  ♦  c 

C*k*7xv»7^,  (12) 

♦  2 

0  *  7  •  V  »  7  X  ,  (13) 

where  D  Is  the  horizontal  divergence.  In  spherical  coordinates,  we 
can  express  (11)  In  component  form 

u  .  _  1  3*  +  1  ,  (14) 

a  a  cos  $  ax 

v  -  ♦  liH  .  (15) 

a  cos  0  ax  a  a<t> 


where  u  Is  the  zonal  velocity  component;  v  Is  the  meridional  velocity 
component,  a  Is  the  earth's  radius;  ♦  Is  the  latitude;  and  x  is  the 
logltude.  If  we  further  define 


U  =  u  cos  $  ,  (16) 

V  =  v  cos  *  ,  (17) 


It  follows  from  (14)  and  (15)  that 


,j  .  -  cos  »  a»  +  1  ax  .  (18) 

a  a*  a  ax 

y  ■  1  a»  ■>.  cos  »  ax  ,  (19) 

a  ax  a  3$ 

♦ 

We  can  also  express  the  frictional  force  F  In  component  form  as 


F  ■  F  1  +  F  j  (20) 

X  4 


i  * 


where  F  and  F  are  the  longitudinal  and  mendlona1  components, 

'  » 

respectively.  The  parameteri zati on  of  these  components  <iH  pe 
detailed  later  in  this  chapter. 

In  order  to  control  spectral  blocking  ( P  u  r  1  and  bourke, 
that  results  from  horizontal  truncation,  diffusion  is  applied  to  the 
vorticity  equation.  This  technique  dampens  deviations  *ro*  equili¬ 
brium  fields  but  not  the  equilibrium  fields,  a i tn  this  specification 
added  and  using  (12),  (14)  to  (20,  we  can  express  the  vorticity 
equation  for  the  troposphere  as  follows: 


2 

37  V  a 

at 


I  [  £*  ♦  COS  *  1®.  ] 

r  n  a> 


2 

K  7  ( 

h 


7  I  -  . 


e  i 


a  cos  « 


'  c .  i 


where  we  define 


A  s  (c  ♦  flu  ♦  o ( 3 v/5o 1  ♦  ( RT  1  /a ) cos  * { 5 q / 3 ♦ )  -  r  cos  t  (2? 

* 


8  :  (t  *  f)v  -  a( 3u/3o 1  -  (RT ’ /a) ! 3q/a a)  ♦  F  cos  t 

x 


22) 


where 


q  =  ln(o  P  ♦  P  ) 
SR 


Is  the  horizontal  diffusion  coefficient.  The  e  subscript  represents 
equilibrium  vorticity  at  the  start  of  model  integrations. 

For  the  upper  atmosphere,  the  parallel  vorticity  equation  may  be 
expressed  as  follows: 

2  A  ^  2  2  ( r  j 

iZ_i  •  -  1  ,  [  ii  +  cos  *  il  ]  ♦  K.  [  7  (;-c  )  ♦  ]  (24) 

at  5T  ax  34  h  e  2 


a  cos  4 


a 
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where 


A  *  ( ;  »  f)U  ♦  w(3v/'-z)  -  F  COS  $  ,  (25) 


3  *  (;  ♦  f)V  -  w(3u/3z)  ♦  F  cos  $  ,  (26) 

A 


for  w  r  32  dt. 

lombining  (21)  and  (24)  into  a  single  expression,  we  can  write 


i;  *t  *  2  ♦  K  [t2 (c-c  )  *  2(;-c  )/a2]  (27) 

n  e  e 

wnere  2  takes  on  the  values  of  the  non-dl ffusi ve  terms  on  the  right 
hai<j  side  of  (21)  for  the  troposphere  and  the  right  hand  side  of  (24) 
for  the  stratosphere  and  mesosphere. 

It  can  be  shown  by  taking  horizontal  divergence  of  (1)  that  the 
divergence  tendency  equation  takes  the  following  form: 

2112*1  »  7  •  [-  (;  ♦  f)I  X  v]  -  72(t  +  »  RT) 

n  q 

-7-'^T'7q-*-3(3v/3o)-F].  (28) 


Expanding  (29)  in  spherical  coordinates,  introducing  horizontal  dif¬ 
fusion,  and  grouping  terms  as  with  the  vorticity  equation,  we  obtain 

2  *  »  f  — 

V  ‘-St  -  cos$  -  v  (E  ♦  *  +  RTq) 
a  cos  * 


r  2  2 
♦  <  [7  D  ♦  2 ( D/a  )] 
h 

where 

2  2 

C  ;  V  •  V  ,  U  ♦  V 

w  -  —5—  - r 


(29) 

(30) 


2  cos  t 


For  the  upper  portion  of  the  model,  we  have 

TT*  ■  — c  It  -  cos*  It  3  •  ,?IE  *  •’ 
a  cos  $ 

2  2 

+  K  [7  D  ♦  2 ( D/a  )]  .  (31) 

h 

As  a  single  combined  equation,  we  can  write  the  divergence  tendency 
equation  as 

3D/3t  -  P  -  72{a  ♦  cRTq)  +  k  [v2D  ♦  2(D/a2)]  ,  (32) 

h 

where  c  *  0  for  the  upper  portion  of  the  model  and  c  ■  1  for  the 
troposphere  and  where  P  takes  on  the  remaining  values  on  the  right 
of  (29)  that  are  not  explicitly  shown  In  (32)  and  on  the  right  of 
(31)  In  a  similar  fashion  for  the  lower  and  upper  portions  of  the 
model  respectively. 

Expanding  the  thermodynamic  equation  (2)  and  expressing  tem¬ 
perature  In  terms  of  a  layer  mean  and  deviation  al’ow  us  to  write 
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i  r»ur  ♦  cos  *  ^T'i  ♦  or  ♦  a uht-.t*) 

2  ax  a#  • 
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h  «  V  T 
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[  f  (p  o+v*vp  )do 
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aV  *7P  ♦ 


m  p 
R  R  ] 
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p  *  a 
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*  P. 


RT  [  f  ( V *7P  )do  -  aV*7P  ♦  P  5’  R.RJ 

w  *  f  #  #  j 

c  (oP  *p  )  0  " 


"  p. 


RT  fl* [  /  Ode  *  7  R  R  1 


Cp  0 


la  _  -  TOP 

-  (3T/ao)[a/  Oda  -  J  Dda  ♦  7(a-l  )2L  exp(JL)  /  0  exp(i£)d2] 

H  Z 


o  n  2  T0P 

.  RT  J  [  /  Oda  ♦  7  Hi  expGl)  /  0  exp(i£)dz]  ,  (33) 

CO  H  H  Z  R 

p  R 

and  the  thermodynamic  equation  for  the  upper  atmosphere 


ar  .  -  l  [  aum  +  cos  *  3 < yT '  Vl  ♦  DT '  -  W  *n 


7T 


7T 


IT 


a  cos  7 
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where  a(z)  Is  the  Newtonian  heating-cooling  coefficient  given  by 

Holton  (1976),  yfj  denotes  the  vertical  diffusion  term, 

s  i  1  *  6'  •  ( 1  /P  ),  i  ;  Til  ♦  S'  *  P  A  <jP  ),  the  overbar  indl- 
*  •  •  ft 

cetes  the  horizontal  domain  average,  cp  is  the  specific  Heat  of  dry 
air  at  constant  pressure,  and  N?  Is  the  buoyancy  f requency  assumed  to 
be  constant  throughout  the  stratosphere  and  mesosphere.  The  term 
RT'w/CpH)  In  (34)  was  neglected  In  the  foraulatlon  of  Lordl  et  al . 
(1980). 

The  surface  pressure  tendency  Is  determined  froa  the  tropospheric 
continuity  equation  In  the  fora 

1  ♦  w  p 

ftq  .  -o  f  (  (V  •  9P  Ido  ♦  P  i  *  *  1 

IT  eP"”p "  J  *  * 

•  R  ° 


1  w  p 

-  06' [  l  Ode  *  T  *  B  ] 

0  H 

1  D„  Z  Zj 

-  a!  [  J  Odo  ♦  T  JL  exp(JL)  j  0  e*p(.I*JdZj  .  (35) 

0  H  HZ  H 

ft 

where  q  »  tnloP  +p  ).  In  Eqs.  (33)  -  (35),  we  have  Isolated  dlver- 
*  ft 

gence  with  layer  aean  terms.  Indicated  by  the  overbar,  and  tlae 
varying  terms.  Indicated  by  the  prime,  so  that  we  will  be  able  to  use 
seal-implicit  time  differencing. 

The  hydrostatic  equation  can  be  written  In  terms  of  layer  aean 
quantities  and  deviations  froa  layer  means.  For  the  troposphere,  we 
have 


34/ 3q  -  -  RT, 


(36) 
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3*‘/5q  -  -  RT'  .  (37) 

For  the  upper  atmosphere  In  terms  of  deviations  from  the  layer  means, 
we  can  write 

ae'/3z  -  RT'/H  .  (38) 

For  the  upper  part  of  the  model,  the  frictional  components  F^ 

and  F  are  represented  by  the  Rayleigh  friction  parameterization. 

In  the  tropospheric  part  of  the  model,  vertical  diffusion  Is  also 

parameterized  In  these  components  along  with  Rayleigh  friction. 

Following  the  vertical  formulations  of  Bourke  et  al .  (1977),  we 

♦ 

can  express  the  tropospheric  components  of  F  as  follows: 


F  -  (flj  !Li  -  F  (U  -  u  )  , 

\  p  TcT  r  e 

(39) 

F,  V  • 

(40) 

where  F  is  the  Rayleigh  friction  coefficient,  u 
r  e 

equilibrium  zonal  wind  components;  and 

Is  the  Initial 

V  75  • 

(41) 

2 

T  •  p  (S_)K  *1  » 

a  p  v  3o 

(42) 

for  (tensity  p  and  vertical  diffusion  coefficient  Kw  defined  as 

K  '  oil.)  v2|”J 
v  p  '5a' 
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Y 


m 

t 


m  1mA 

P  (u)« 

t 


(50) 


u  »  sin  ♦  .  (51) 

m 

P  (u)  Is  an  associated  Legendre  polynomial  normalized  to  unity  with 
m,  the  longitudinal  wavenumber,  and  t,  the  latitudinal  Index.  L 
determines  the  limit  of  the  paral lelogramlc  truncation.  The  {} 
terms  on  the  right  of  (46)  -  (49)  represent  the  respective  spherical 
harmonic  coefficients.  J  Is  the  longitudinal  wavenumber  truncation. 

After  the  vertical  derivatives  and  non-linear  products  have  been 
formed  at  grid  points,  the  resulting  terms  can  be  transformed  In 
terms  of  Fourier  series  as  follows: 

1mA 

(2.  0.  T,  P,  Y'[0],  r[T‘ ]}  «  t  {Z  ,  p  ,  T  ,  P  ,  9  ,  <  }e  (52) 

IJjSaJ  Ifl  H  W  IW  W 

where  the  m  subscript  represents  the  respective  Fourier  coefficient. 
The  Legendre  transform  defined  by 

m  m 

(5  ■  /  (  )  P  (u)cosa  da  (53) 

t  -w/2  "»  t 

can  be  applied  to  the  Fourier  coefficients  In  (52)  to  obtain  the 
appropriate  spherical  harmonic  coefficients.  Since  In  the  model, 
the  Fourier  series  defined  by  (53)  are  formed  at  Gaussian  latitudes, 
(53)  can  be  computed  exactly  for  each  term  up  to  the  point  of  trunca¬ 
tion. 


Initial  fields  of  geopotentlal ,  temperature,  stream  function  and 
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log-pressure  must  be  specified  before  the  start  of  model  Integrations. 
For  these  Initial  conditions,  we  specify  zonal  profiles  that  are  non- 
dlvergent  and  essentially  have  a  non-dl vergent  tendency.  These 
latter  conditions  are  essential  so  that  high  frequency  oscillations, 
which  could  adversely  affect  the  model's  behavior,  are  basically 
eliminated. 

To  arrive  at  a  balanced  state  In  the  troposphere,  we  Initially 
started  with  constant  pressure  level  zonal  temperature  deviations 
from  level  means,  based  on  January  climatic  tables  from  Oort  and 
Rasmusson  (1972).  Data  for  polar  latitudes  were  extrapolated  from 
January  1974  and  1976  data  from  the  National  Meteorological  Center 
(NMC)  observational  grids.  The  above  temperature  fields,  which  were 
specified  at  5°  latitude  Increments,  were  first  linearly  Interpo¬ 
lated  to  Gaussian  latitudes  used  In  the  model. 

Averaging  the  diagnostic  level  geopotentials  to  obtain  prognos¬ 
tic  level  values  as  Is  done  In  the  model  and  using  the  corresponding 
temperature  and  pressure  fields,  we  can  then  quadratlcally  solve  the 
gradient  wind  balance  equation  for  our  o-coordinate  system 

u(f  +  u  tan$)  «  _  1_  .  qR(T  +  T1 )  1_  (54) 

a  a  3$  a  3$ 

*  R 

for  the  zonal  velocity  field  u.  The  mean  temperature  profile  T  Is 
assumed  to  be  that  of  the  standard  atmosphere  at  45°  N. 

The  orographic  forcing  In  the  model  Is  very  Idealistic  and 
simplistic.  The  topographic  pattern  was  placed  In  such  a  manner  as 
to  simulate  the  continental  land  mass  distribution  of  the  Northern 


Hemisphere  (Bourke).  The  longitudinal  distribution  of  the  surface 
height,  z^,  at  45N  can  be  expressed  empirically  as  follows: 
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zJ45N)  «  A 


s1n(m  X)  +  A  s1n(m  x  -  *  )  +  z  , 
1  2  2  2  * 


(55) 


where  A  and  A  are  the  amplitudes  of  zonal  wave  numbers  m  •  1  and 
12  1 

m  >2,  respectively;  7  Is  the  Initial  flat  surface  height  and  was 
2  * 

used  to  determine  the  Initial  mean  surface  pressure  p^  based  on  stan¬ 
dard  atmosphere  profiles.  The  latitudinal  variation  can  be  expressed 
by 

2 

zje)  •  zJ45N)  sin  (2$)  .  (56) 


The  model  which  generated  the  data  base  for  this  thesis  had  a 

7  *  A  ■  240  m  and  A  *  150  m.  The  maximum  z  for  the  simplistic 
*  2  1  * 

European/Aslan  land  mass  Is  630  m  (A  ♦  A  +  7  )  and  for  the  North 

1  2  * 

American  continent,  z  was  330  m  (A  -  A  +  7  ). 

*  2  1  * 

To  prevent  the  model  from  being  shocked  by  the  sudden  Intro¬ 
duction  of  topography,  It  was  slowly  turned  on  according  to  the 
following  algorithm: 


♦'  *  g(z  -  z  )[1  -  exp(-t/t  )]  ,  (57) 

*  *  *  o 

5 

for  time  t  where  t  *  2.5  x  10  s.  With  this  formulation,  ♦  attains 
o  * 

about  90%  of  Its  maximum  value  after  7  days. 

We  have  tested  the  spectral  model  and  Integrated  for  forty  days 
to  obtain  stable  solution.  Preliminary  reports  obtained  from  the 
numerical  simulation  with  the  use  of  this  model  compare  favorably 
with  the  observations. 
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CHAPTER  3 


ANALYSIS  OF  THE  BEGINNING  FORCES  ON  THE 
DEVELOPMENT  OF  THE  STRATOSPHERIC  WARMING 

The  31  layer  model  was  Initialized  with  mean  zonal  winds  as 
depicted  In  Fig.  2a  for  case  1  (Cl)  and  Fig.  2b  for  case  2  (C2),  and 
with  orographic  forcing  as  shown  In  Fig.  3.  For  nearly  all  quanti¬ 
ties  computed  In  the  analysis  of  this  model,  C2  will  have  magnitudes 
two  to  six  times  as  large  as  Cl.  Since  the  winter  climatology  Is 
adjusted  to  give  a  non-dlvergent  flow  and  essentially  a  non-dl vergent 
tendency,  there  are  no  perturbations  developing  until  the  orographic 
features  are  Introduced.  This  means  that  the  Initial  perturbations 
of  the  model  will  be  those  caused  by  the  orography  and  thus  kfV  Is 
the  Initial  perturbations  of  the  model.  The  Induced  orography 
(Fig.  3)  Is  fairly  simple  with  Its  maximum  height  at  45°N.  Yet,  the 
maximum  flux  of  ffTTr  after  day  6  occurs  near  70°N. 

In  referencing  various  levels  of  the  atmosphere,  we  will  use 
approximately  13.5  Km  as  the  tropopause,  47  Km  as  the  stratopause  and 
85  Km  as  the  mesopause.  These  values  are  based  on  the  vertical 
temperature  profile  for  the  U.S.  Standard  Atmosphere  (1976).  For  Cl 
we  selected  day  18  and  day  30  as  representative  charts  of  the  zonal 
mean  latitude-height  distributions.  Day  18  represents  the  general 
configuration  of  the  quantities  before  the  wind  reversal  (major 
warming)  and  day  30  the  configuration  after  the  warning.  These  days 


are  also  the  period  of  wave  transition  In  the  lower  levels  of  the 
model.  For  C2  we  selected  day  26  and  day  30.  Day  26  Is  representa¬ 
tive  of  conditions  prior  to  and  during  the  first  warming  period, 
while  day  30  represents  some  of  the  significant  changes  that  occur 
during  the  last  stages  and  just  after  a  warming  period.  We  will  also 
comment  on  day  34  of  C2  since  It  seems  to  reflect  a  transitory  condi¬ 
tion  between  warming  surges. 

3.1  Effect  of  the  Surface  Forcing  Function 
on  the  Development  of  Waves 

By  day  2,  after  the  Introduction  of  the  forcing  function,  there 
Is  an  accumulation  of  horizontal,  as  well  as  vertical  velocity 
variances  at  the  top  of  the  model.  However,  by  day  6  the  variance 
has  shifted  from  the  top  of  the  model  to  the  stratosphere.  After 
this  shift  of  variance,  It  appears  that  most  of  the  development 
occurs  In  the  region  of  the  stratosphere.  By  looking  at  the  zonal 
mean  vertical  propagation  of  geopotential  W and  the  velocity 
variances  that  it  Induces,  we  try  to  develop  a  better  understanding 
of  the  mechanisms  Involved  that  bring  about  the  rapid  stratospheric 
warmings. 

On  day  18  of  Cl,  there  is  a  maximum  of  W at  about  76°N  and 
50  Km  (Fig.  4a).  By  day  30,  we  note  the  flux  of  geopotential  has 
dropped  to  about  one  third  the  magnitude  of  day  18  and  moved  slightly 
south  indicating  a  more  quasi -stationary  state.  On  day  26  of  C2  a 
similar  location  of  W ' * '  to  Cl's  day  19  (Fig.  4c).  However,  by  day  30 
of  C2,  a  cell  of  sinking  geopotential  has  develooed  In  the  area  for 
formally  rising  geopotentlal  on  day  26  (Fig.  4d),  but  by  day  34  the 
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model  has  returned  to  a  near  pre-warming  configuration.  The  differ¬ 
ence  of  magnitudes  from  day  18  to  day  30  In  Cl  and  from  26  to  day  30 
of  C2  Indicates  that  the  vertical  energy  propagation  has  been  reduced. 
Since  our  forcing  function  Is  constant,  this  Implies  that  there  has 
been  a  dynamic  change  In  the  environment.  The  contrast,  however,  In 
the  difference  of  relative  magnitudes  on  day  30  to  day  18  of  Cl  and 
day  34  to  day  26  of  C2  Is  our  first  Indication  of  a  major  difference 
In  the  development  of  the  two  cases.  It  seems  that  C2  and  Cl  are 
undergoing  slightly  different  atmospheric  changes.  These  differences 
of  dynamic  effects  seem  to  be  primarily  In  the  area  of  wave  to  wave 
and  wave  to  mean  Interactions.  In  an  effort  to  discern  these  critical 
differences  a  limited  spectral  analysis  of  the  two  cases  was  accom¬ 
plished.  The  results  of  this  analysis  is  discussed  in  Chapter  6. 

The  consequences  of  the  vertical  flux  of  geopotential  can  be 
seen  most  readily  In  terms  of  variance.  As  noted  earlier,  the 
concentration  of  the  variance  of  the  velocities  began  at  the  top  of 
the  model,  then  shifted  to  the  stratosphere,  primarily  In  the  polar 
regions,  by  day  6  In  both  Cl  and  C2. 

It  should  be  noted,  however,  that  u'z  had  two  additional  maxlmums, 
besides  the  one  In  the  polar  stratosphere.  One  of  these  was  at  mid¬ 
latitude  In  the  stratosphere  and  the  second  at  25°N  at  the  tropopause. 
This  last  variance  of  u  Is  the  dominant  maximum  maximum  by  nearly 
twice  the  stratospheric  quantities.  The  vertical  velocity,  w,  had  an 
area  of  significant  variance  located  at  the  tropopause  and  about  35°N. 
C2  Is  very  similar  to  Cl  during  the  first  ten  days  of  development 
except  some  of  Its  maxima  of  variances  are  shifted  northward  about  5 
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degrees . 

As  the  model  progresses  with  time,  we  can  see  from  Fig.  5  and 
Fig.  7  that  the  variances  are  concentrated  at  the  stratopause.  By 
comparing  Fig.  5  and  Fig.  6,  It  Is  fairly  easy  to  see  that  the  maximum 
variance  of  7  occurs  In  the  areas  of  maximum  gradient  of  u”.  Thus  for 
Cl  there  are  two  maxlmums  of  variance  on  day  18  (Fig.  6a)  with  the 
northern  maximum  dlmlnlsning  by  day  30  as  the  polar  night  jet  (PNJ) 
decreases  In  velocity  and  reverses.  In  C2,  however,  we  note  that  the 
symmetry  of  variance  magnitudes  (Fig.  5c)  on  day  26  shifts  to  a 
dominant  cell  at  the  polar  region  by  day  30  even  though  the  PNJ  has 
diminished.  This  could  be  due  to  the  model  cycling  back  Into  a  single 
core  jet  by  day  34. 

Relating  the  variance  of  IT  to  that  of  7,  the  meridional  velocity, 
there  are  three  features  which  stand  out.  First,  there  tends  to  be 
only  one  area  of  maximum  of  v'2  as  opposed  to  lT7^  having  two.  The 
second  Is  that  v'^  is  between  the  areas  of  maximum  u'2  and  the  third 
Item  worth  noting  Is  that  u7^  and  T7^  have  a  very  large  maximum  on 
day  30  In  C2,  whereas,  In  Cl  u'2  peaks  on  day  30,  but  v'^  Is  decreas¬ 
ing  on  day  30. 


3.2  Comparing  the  Vertical  Flux  of  Geopotential 
to  the  Meridional  Flux  of  Zonal  Momentum  and 
the  Variance  of  the  Zonal  and  Meridional  Velocity 

Looking  back  at  day  18  and  day  30  for  Cl,  it  appears  that  the 


maximum  component  of  vertical  flux  of  geopotential  occurs  at  approxi¬ 
mately  the  same  height,  40  km,  as  the  maximum  of  the  meridional 
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flux  of  turbulent  momentum  ( u ' v ' ) .  This  Is  not  too  surprising 
since  this  Is  also  the  approximate  height  of  the  variance  maxima  of 
u"  and  7.  Comparing  Fig.  8  to  Figs.  5  and  7,  we  also  note  that  the 
maximum  of  u’v'  Is  between  the  maximum  centers  of  u'Z  and  v'2  except 
on  day  30  of  C2  where  the  maximum  of  u'2  and  v’2  are  both  near  90#N. 

The  area  of  maximum  u'v'  In  Cl  Is  centered  at  55°N  on  day  18  and 
45°N  by  day  30. 

C2  Is  very  similar  on  day  26  to  day  18  of  Cl,  but  day  30  of  C2 
Is  quite  different.  The  area  of  strong  u'v’  convergence  on  day  26 
has  been  replaced  with  a  much  weaker  area  of  convergence  and  an  area 
of  divergence  has  developed  at  72°N  and  43  Km.  This  seems  to  be  a 
transient  phenomenon  which  disappears  by  day  34.  Many  of  these 
momentum  transport  changes  In  Cl  and  C2  seem  to  be  related  to  phase 
changes  of  the  spectral  components.  These  changes  will  therefore  be 
discussed  In  more  detail  In  the  chapter  on  spectral  analysis. 

3.3  Relationship  of  the  Vertical  Flux 
of  Geopotential  to  the  Meridional  Flux 
of  Sensible  Heat 


The  linearized  relationship  [Eq.  (58)]  developed  by  Charney  and 
Drazln  (1961)  and  Ellasen  and  Palm  (1961)  relates  vertical  transport 
of  geopotential  to  meridional  transport  of  sensible  heat  and  can  be 
viewed  as  Eq.  (59). 


f 

r  (U-C) 


3^ 

3 X  3Z 


(58) 


wT  *  W  (U-C)  . 


(59) 
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This  relation  Indicates  that  the  larger  the  amplitude  of  wTTr  the 
larger  the  amplitude  of  v’T' .  However,  as  we  compare  this  linearized 
equation  to  the  spectral  model  diagrams,  we  find  only  a  degree  coin¬ 
cidence.  For  example,  day  18  of  Cl  (Fig.  9a)  shows  that  v’T^x  Is 
south  of  w^max  (^9*  4a)-  This  Is  also  tve  case  on  day  26  of  C2. 

By  the  end  of  the  warming  period,  the  meridional  transport  of 
sensible  heat  and  vertical  transport  of  geopotentlal  come  close  to 
coinciding.  This  Indicates  that  the  linearized  relationship  does  not 
hold  well  for  highly  nonlinear  (nonstationary)  state,  but  may  apply 
for  a  quasi -stationary  condition.  Another  feature  which  casts 
questions  on  the  applicability  of  the  linearized  equation  Is  that 
within  the  area  of  positive  v'T’  flux  on  day  30  of  C2  (Fig.  9d)  there 
Is  an  area  of  -  w’*'  (Fig.  4d).  The  only  way  that  this  would  be 
possible  would  be  that  C  exceeds  u.  On  day  30  of  C2  u"  has  a  value  of 
approximately  55  m/sec  at  this  location.  It  Is  unlikely  that  C  would 
exceed  this  value. 
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CHAPTER  4 


ANALYSIS  OF  THE  ZONAL  MEAN  CIRCULATION 
IN  THE  UPPER  ATMOSPHERE 

Recall  that  even  though  the  highest  point  of  the  orographic 
forcing  was  at  459N,  the  strongest  flux  of  w' T*  in  the  stratosphere 
was  at  about  70°N  on  day  18  of  Cl  and  day  26  of  C2.  The  latitude- 
height  chart  of  w" t'  on  day  3  shows  that  the  maximum  flux  of  w'$ '  is 
near  35°N  in  the  troposphere.  As  this  eddy  flux  propogates  Into  the 
stratosphere,  it  shifts  northward  to  about  65°N.  As  the  vertical 
flux  of  geopotential  reaches  a  maximum  in  the  stratosphere,  the  mean 
zonal  meridional  circulation  also  peaks. 

4.1  Mean  Meridional  Circulation 

Examining  Figs.  10a  and  10b  for  Cl,  we  can  see  the  zonal  mean 
meridional  circulation  as  noted  by  Matsuno  (1971).  This  circulation 
was  established  by  day  10  of  the  model  Integration.  The  circulation 
consists  of  rising  motion  in  the  polar  region  across  the  stratopause, 
southward  motion  in  the  lower  mesosphere  and  sinking  motion  into  the 
stratosphere  at  mid-latitudes.  This  circulation  coincides  with  the 
area  of  maximum  fluxes  of  geopotential,  sensible  heat,  and  momentum. 

An  interesting  feature  of  this  circulation  is  that  it  extends 


from  the  mesosphere  to  the  troposphere  such  that  the  northerly  portion 
of  the  cell's  circulation  is  confined  primarily  to  the  troposphere. 
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From  day  13  to  day  30  of  Cl  there  has  been  a  weakening  of  the  zonal 
mean  meridional  circulation  and  Its  center  has  shifted  south  from 
65°N  to  55°N.  With  Figs.  10c  and  lOd  the  evolution  of  the  mean 
meridional  velocity  7max  can  be  followed.  Figure  10c  Is  the  maximum 
of  v  plotted  every  fourth  day  to  show  changing  mean  amplitudes  with 
time.  Figure  lOd  Is  a  latitude-time  section  taken  at  40.5  km  which 
depicts  7  fluxuations  In  both  latitude  and  amplitude.  The  first 
maximum  of  v  occurs  on  day  14  and  the  second  on  day  26  which  holds 
through  day  30  with  a  minimum  on  day  22.  It  should  also  be  pointed 
out  that  by  day  34  the  circulation  cell  Is  virtually  dissolved. 

To  complete  the  meridional  circulation  analysis  shown  In  Figs. 
10a  and  10b,  we  will  look  at  the  evolution  of  the  vertical  velocities 
with  time.  Figure  11a  shows  wmax  (solid  curve!  and  wm^n  (dashed 
curve)  respectively  of  the  northern  and  southern  branch  of  the  zonal 
mean  meridional  circulation  and  Fig.  lib  depicts  w  as  a  function  of 
time  similar  to  Figs.  10c  and  lOd  respectively. 

It  Is  obvious  that  Fig.  11a  has  a  geometric  configuration  simi¬ 
lar  to  Fig.  10c  except  that  Fig.  11a  has  the  mirror  effect  of  the 
positive  and  negative  vertical  velocity  maxima.  These  fluxatlng 
maxima  Indicate  a  pulsing  meridional  circulation. 

For  C2  we  find  that  the  mean  meridional  circulation  on  day  26  is 
very  similar  to  Cl  (Fig.  12a),  but  day  30  has  a  second  circulation 
consisting  of  a  northward  maximum  of  7  at  75°N  and  40  Km  and  sinking 
motion  near  90°N.  This  cell  is  of  very  short  duration  and  is  not 
depicted  on  the  7max  evolution  charts  (Figs.  12c  and  12d)  or  the  w 
evolution  charts  (Fig.  13a  or  13b).  The  C2  evolution  charts  of  the 
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mean  meridional  maxima  show  three  distinct  maximum  values  of  v.  The 
first  two  surges  of  meridional  circulation  occur  on  the  same  day  as 
Cl.  However,  Instead  of  the  meridional  circulation  going  to  near 
zero  on  day  34  as  In  Cl,  It  again  Increases  through  day  40.  Because 
of  this  Increase  In  activity  on  day  40,  the  model  was  Integrated  by 
Koermer  through  day  50.  This  Integration  showed  that  day  40  was  the 
last  day  of  significant  Increase.  A  second  difference  between  Cl  and 
C2  Is  that  In  Cl  the  meridional  circulation  began  to  shift  south 
after  day  26,  whereas,  In  C2  the  meridional  circulation  continued 
moving  northward  through  day  40  (Fig.  12d  and  Fig.  13b).  In  Fig.  13a 
we  again  note  a  similar  mirror-type  Image  of  the  vertical  velocity 
magnitudes.  C2  also  shows  an  approximate  8-12  day  cycle. 

4.2  Theorems  of  Noninteraction 

Figure  14  is  a  graphic  display  of  the  relationships  between  the 
zonal  mean  v'T'  (solid  curves)  and  w  (dashed  curves).  The  relations 
given  by  Eqs.  (53)  and  (54)  represent  time  averaged  relations  (3) 
between  the  zonal  mean  meridional  circulation  and  the  flux  of  sensible 
heat  and  the  transport  of  zonal  momentum. 

a  •  ^  a 

w  *  _J _ [  -  1  _ .  JL  (vT  cos  $)+£_]  (53) 

*  a  cos  4  3<$>  u 

Y  -Y  P 
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a  •  —  a  a  a 

v  *  4  [  1  !_  (vu  cos  <f> )  +  i_(wu)  -  tan  »fvu)  _  f  ] 

r  a  cos  <j>  3p  a  1 

(54) 

From  Eq.  (53)  we  would  expect  rising  motion  to  be  In  an  area  of 


horizontal  convergence  of  v'T' .  Likewise,  In  an  area  of  horizontal 
divergence  of  v'T'  there  should  be  sinking  motion.  Figure  14a  Is  a 
representative  height-latitude  chart  depicting  this  relationship. 

The  relationship  seems  to  hold  true  for  both  Cl  and  C2.  Equation 
(54)  shows  the  relation  between  v  and  the  momentum  fluxes  of  uTv'r  and 
w'u'.  In  the  relation  to  v'u' ,  v  should  be  negative  In  the  area  of 

3 

negative  gradient,  as  well  as  the  maximum  negative  portion  of  v  at 

the  center  of  v’u' .  Figure  14b  Indicates  that  1  (v'u'  cos  ♦)  1 

'  _  a  cos  » 

the  dominant  of  the  v'u’  terms.  The  w'u  *  term  of  Eq.  (61)  should 
give  meridional  flow  In  the  area  of  maximum  vertical  gradient. 

Figure  14c  Indicates  that  while  w’u'  Is  a  possible  factor  It  is  not 

3 

the  dominant  Influence  on  v. 

It  may  be  pointed  out  from  the  following  equation  (Kao,  1981). 

wT  *  CpT  |  _  1  3_  (7r  cos  »)  -  !_  (wf)  1  (62) 

g  a  cos  *  a*  az 

and  Eq.  (53)  that  the  region  of  zonal  mean  vertical  flux  of  tempera¬ 
ture  generally  coincides  with  that  of  zonal  mean  vertical  velocity. 


CHAPTER  5 


EVOLUTION  OF  THE  ZONAL  MEAN  MAXIMUM  OF  ZONAL 
VELOCITY,  TEMPERATURE,  AND  CROSS  CORRELATIONS 
OF  u,v,w,  and  T. 

Since  the  difference  between  a  major  and  a  minor  warming  In  the 
stratosphere  depends  on  the  characteristic  of  polar  night  jet  (PNJ), 
we  will  concentrate  on  the  evolution  of  this  jet.  According  to  the  WHO 
definition.  If  the  PNJ  mean  winds  reverse  from  westerlies  to  easterlies, 
then  It  Is  considered  a  major  warming.  All  other  strong  warmings  are 
classified  as  minor  warmings. 

5.1  Evolution  of  Zonal  Mean  Winds 
and  Temperatures 

In  both  cases  we  find  that  the  PNJ  becomes  a  separate  entity  on 
day  14  which  corresponds  to  the  first  surge  of  meridional  circulation. 

For  Cl  (Fig.  15a)  we  note  that  the  PNJ  (solid  line)  holds  steady 
from  day  14  through  day  22,  then  de-accelerates  rapidly  with  the  mean 
winds  reversing  on  day  26.  During  the  period  of  simulation,  the  meso¬ 
spheric  jet  (dashed  line)  has  a  fairly  steady  decline  through  day  30, 
then  gradually  Increases.  The  tropospheric  jet  (dot  line)  and  the 
mid-latitude  easterlies  (dot  dashed  line)  maintain  a  relatively  even 
magnitude  during  the  entire  Integration  period.  See  Fig.  6d  for  jet 
locations. 


For  C2  (Fig.  15b)  the  Intensity  of  PNJ  Increases  sllghty  from  day 
14  through  day  22,  then  de-accelerates.  From  day  30  to  day  34,  the  de- 
acceleration  rate  slows  during  a  transitory  period.  The  C2  mesospheric 
jet  has  a  much  larger  fluxatlon  than  the  Cl  jet.  This  Is  the  opposite 
of  the  findings  of  Koermer  and  Kao  (1980),  In  their  analysis  of  the 
1977  major  and  the  1976  minor  warming.  As  In  Cl,  the  tropospheric  and 
mid-latitude  easterlies  maintain  relatively  steady  velocities. 

The  evolution  of  the  changes  of  the  zonal  mean  temperatures  maxima 
of  Cl  are  depicted  In  Fig.  16.  Figure  16a  for  Cl  shows  a  nearly  linear 
Increase  In  the  stratospheric  temperature  maximum  (solid  curve)  from 
day  10  through  day  34,  whereas  the  mesospheric  temperature  maximum 
experiences  a  steady  decline  through  the  simulation  period.  This 
indirectly  shows  that  there  Is  some  law  of  compensation  between  the 
temperature  changes  In  the  stratosphere  and  those  In  the  mesosphere. 

For  C2  we  have  a  slightly  different  warming  transition  In  the 
stratosphere,  with  a  cooling  period  from  day  30  to  34,  and  the  final 
warming  on  day  38.  To  a  large  degree,  we  see  a  mirror  effect  of  this 
stratospheric  warming  as  cooling  In  the  mesosphere. 

The  height  of  the  warming  maximum  for  Cl  Is  fairly  steady  (solid 
line)  (Fig.  16b),  with  some  sinking  during  about  the  last  half  of  the 
period.  C2  (Fig.  17b),  shows  more  oscillation  in  the  vertical  position 
of  maximum  warming  (solid  line)  with  Its  position  rising  through  day 
30,  then  sinking.  The  height  of  the  maximum  cooling  (dashed  line) 
changes  quite  dramatically  (Fig.  16b),  from  day  14  *•-  day  18  In  Cl  and 
then  remains  quasi -stationary  with  some  shift  In  height  simultaneously 
with  the  stratosphere  warming  height.  C2  cooling  height  (dashed  line) 
(Fig.  17b)  does  not  have  the  rapid  height  change  of  Cl,  but  rises  dally 
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from  day  14  through  day  34.  The  sinking  of  the  cooling  maximum  In  the 
mesosphere  and  the  downward  shift  of  the  stratospheric  warming  Is  not 
as  simultaneous  In  C2  as  In  Cl.  The  latitude  of  the  stratospheric 
maximum  warming  position  Is  closely  correlated  to  the  periods  of  maximum 
warming,  especially  In  C2.  These  maximum  temperature  change  areas  are 
also  closely  related  to  the  de-acceleration  of  the  PMJ  winds.  As  the 
area  of  maximum  temperature  change  shifts  north,  the  magnitude  of  the 
warming  increases.  There  also  seems  to  be  some  correlation  between 
the  rate  of  movement  north  and  the  rate  of  temperature  change  In  the 
stratosphere.  This  holds  true  for  both  Cl  and  C2.  In  C2,  as  the  area 
of  maximum  warming  moves  south  from  day  26  to  day  34,  It  cools  nearly 
18"C.  However,  as  It  shifts  northward,  day  34  to  day  38,  It  warms 
again. 

The  mesosphere  cooling  maximum  shifts  northward  from  day  10 
through  day  18  for  Cl  (Fig.  16c),  and  exhibits  very  little  latitude 
changes  thereafter.  The  C2  mesosphere  cooling  maximum  has  character¬ 
istics  similar  to  Cl  except  It  shifts  northward  from  day  22  to  day  26. 

5.2  Evolution  of  Fluxes  from  Magnitude 
Time  Charts  and  the  Effect  of  Meridional 
Fluxes  on  Zonal  Mean  Temperature  and  Hind 

Since  the  vertical  flux  of  geopotential  Is  the  first  perturbation 
to  be  formed  from  the  orographic  forcing,  we  will  look  at  the  evolution 
of  the  maximum  values  of  this  term  first. 

In  Cl  the  vertical  flux  of  geopotentlal  has  Its  maximum  values 
In  the  mid-latitudes  shifting  northward  with  time  (Fig.  18b).  When  the 
amplitude  of  7T  reaches  Its  peak  on  day  18  (Fig.  18a),  the  latitude 
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position  seems  to  remain  constant  until  the  warming  trend  develops  on 
day  22.  Coinciding  with  the  PNJ  de-acceleration  Is  the  formation  of  an 
area  of  negative  wTTr  (shown  by  an  H).  As  the  warming  progresses,  the 
area  maximum  shifts  south.  C2  exhibits  more  and  larger  fluxatlons,  but 
of  shorter  duration.  As  In  Cl,  the  maxima  shift  from  the  mid-latitudes 
to  the  polar  regions  near  the  beginning  of  the  Integration.  As  the 
amplitude  of  the  maxima  decrease  (Fig.  19a),  they  seem  to  move  south 
(Fig.  19b),  and  the  amplitude  of  the  minima  Increases  (N).  These 
minima  seem  to  be  related  to  the  decrease  In  the  meridional  circulation 
and  warming  In  the  stratosphere. 

The  Cl  meridional  transport  of  mean  zonal  momentum  maxima  (Fig. 

20a)  has  a  very  similar  configuration  to  the  vertical  flux  of  geopoten¬ 
tial  maxima,  except  that  It  reaches  Its  peaks  four  days  earlier.  The 

i 

latitude  center  of  these  maxima  does  not  reach  the  polar  regions,  but 
does  show  a  northward  shift  (Fig.  20b)  reaching  Its  northernmost 
position  on  day  18.  As  with  wTTr,  the  divergent  area  of  u' v'  forms  at 
the  beginning  of  the  warming. 

In  C2  the  zonal  mean  meridional  flux  of  zonal  momentum  maxima 
evolution  configuration  (Fig.  21a),  is  very  similar  to  the  maxima 
amplitude  of  u'v'  Increases,  the  latitude  position  (Fig.  21b), 

i 

Indicates  convergence.  As  divergent  areas  form,  the  latitude  position 
of  the  convergence  maxima  shifts  south  and  the  magnitudes  decrease. 

The  mean  zonal  meridional  flux  of  zonal  momentum  should  have  a 
direct  effect  on  the  mean  zonal  winds  as  shown  In  Eq.  (63).  It  can  be 
shown 
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that  at  68.9#N  the  turbulent  transfer  of  zonal  momentum  and  the 
Corlolus  torque  were  the  dominant  terms  of  the  relationship.  However, 
these  terms  were  of  opposite  sign,  and  thus  nearly  negated  the  effect. 

The  small  feature  which  does  have  a  direct  correlation  to  the 
warming  Is  the  divergent  area  of  u' v'  In  the  polar  mid-stratosphere. 

As  the  mean  zonal  meridional  flux  of  zonal  momentum  changes  sign.  It 
would  enhance  the  Coriolis  effect  on  de-acceleratlng  the  PNO  winds. 

The  vertical  transport  of  zonal  momentum  (w'u' )  has  a  significantly 
different  configuration  from  u' v' .  In  Cl  Its  maximum  magnitudes  peak 
on  day  18  and  then  gradually  decline  to  day  26  with  only  a  minor  peak 
on  day  30  (Fig.  22a).  The  C2  configuration  (Fig.  22b)  at  w'u1  has  a 
similar  time  arrangement  of  maxima  and  minima  to  urvr  through  day  30 
after  which  there  Is  no  correlation.  These  configurations  plus  the 
magnitude  of  the  term  suggests  that  w’u'  does  not  have  a  direct 
Influence  on  the  stratospheric  warming. 

The  mean  zonal  meridional  transport  of  sensible  heat  has  nearly 
the  same  configurations  as  7rvr  for  each  respective  case  (Figs.  23a  and 
24a).  The  major  difference  being  that  v'T'  does  not  have  the  negative 
maximum  in  Cl  or  C2  that  was  evident  In  w' *'  and  u' v' .  It  should  be 
noted,  however,  that  the  evolution  of  the  amplitude  (Fig.  23a),  of  the 
maxima  of  v'T’  In  Cl  has  a  very  similar  configuration  to  the  Cl  u ' v ' 
maxima  evolution.  After  day  30  the  area  of  v'T*  has  dissolved  to  the 
point  that  no  one  center  can  be  distinguished.  In  C2  the  latitude 
position  of  V'T'  does  not  Indicate  that  magnitude  relation  that  was 
observed  In  the  previous  flux  terms.  It  does,  however,  Indicate  a 
gradual  poleward  movement  of  warm  air  which  is  consistent  with  the 
findings  of  Quiroz,  et  aj_.  (1975). 
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Equation  (64)  relates  the  meridional  and  vertical  transport  of 
sensible  heat  to  the  changes  of  temperature  with  time,  where  (  )s  Is 
the  transient  term. 

(jL  «  -  1  L_  (7r  cos  *)  -  i_  (w7?7).  ♦  (iTT7)  .  <64) 

Tt  c  a  cos  4  34  s  3p  5  c  p  s 

5  P 

Using  a  simple  differencing  program  to  find  gradients.  It  was  deter¬ 
mined  that  the  meridional  transport  term  Is  approximately  an  order  of 
magnitude  larger  than  vertical  transport  term  and  thus  probably  more 
Important.  This  Is  Indirectly  confirmed  by  Figs.  16  and  17  which  show 
the  correlation  between  meridional  shift  of  the  temperature  change 
maximum  an  the  magnitude  of  that  change. 

The  evolution  similarity  of  magnitudes  and  frequencies  of  v'T' 
and  u V ,  as  depicted  In  Figs.  20a  and  22a  and  In  Figs.  21a  and  23a, 
suggests  that  a  relation  exists  between  these  fluxes.  Equation  (58) 
relates  w  to  v'T1  and  Eq.  (59)  relates  v  to  v V .  In  the  Indirect 

a  a 

meridional  circulation  as  shown  In  Fig.  10a,  we  know  that  w  and  v  are 
related  by  the  continuity  equation.  Thus  an  Indirect  relation  should 
exist  between  v'T'  and  u'v' . 


CHAPTER  6 


ANALYSIS  OF  SPECTRAL  INFLUENCE  ON  THE  CHANGE 
OF  THE  MEAN  ZONAL  WINDS  AND 
MEAN  ZONAL  TEMPERATURE 

The  diagrams  of  this  chapter  are  based  on  a  Fast  Fourier  Transform 
of  u,  v,  w,  T,  and  *  into  their  spectral  components.  From  this  spec¬ 
tral  data  base,  we  computed  the  spectral  zonal  mean  amplitudes  of  the 
velocity  components  denoted  as  A(  )(k)  with  (k)  representing  the 
wavenumber.  The  cospectral  values  are  designated  as  C(  )(  )(k)  and 
represent  the  zonal  mean  of  the  flux  terms.  When  selecting  the  area 
maxima,  we  focused  our  attention  on  the  polar  region  since  this  Is 
the  region  of  the  stratospheric  warming.  The  latitude-height  charts 
were  selected  to  maintain  continuity  with  previous  chapters  and  to 
highlight  critical  periods  In  the  warming  process  that  were  indicated 
on  the  evolution  charts. 

The  magnitude  evolution  graph  periods  are  based  on  the  time 
periods  of  maximum  zonal  velocity  changes  In  the  polar  stratosphere. 
Thus  for  Cl,  we  selected  days  16  through  30  and  for  C2  days  22  through 
40.  Since  wavenumber  1  and  wavenumber  2  (hereafter  referred  to  as  wave 
1  and  wave  2  respectively)  have  the  dominant  wave  amplitudes  of  this 
study  [Koermer  (I9*.0)]f  we  will  refer  to  wavenumber  three  and  wave- 
number  four  only  when  they  have  significant  magnitudes  compared  to 
wavenumbers  one  and  two. 


6.1  Amplitude  and  Evolution  Analysis 
of  Au(k),  Av(k),  and  A^tk) 
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One  of  the  most  striking  features  of  Cl  Is  that  for  AtJ(k)  (Fig. 

25),  Av(k)  (Fig.  26),  and  Aw(k)  (Fig.  27)  wave  1  dominates  the  meso¬ 
sphere  and  wave  2  dominates  the  stratosphere.  This  holds  true  through¬ 
out  the  warming  period  even  though  wave  1  does  seem  to  propagate  down¬ 
ward  Into  the  stratospheric  polar  region  by  the  end  of  the  warming. 

This  Implies  that  wave  1  has  a  greater  vertical  penetration  ability 
than  wave  2. 

Au(l)  (Fig.  25a)  starts  with  three  wave  centers  symmetrically 
spaced  from  the  equator  to  the  North  Pole.  Av(2)  also  has  three 
centers  of  maximum,  one  of  which  is  located  in  the  lower  stratosphere. 
This  suggests  that  wave  2  may  have  more  Interaction  with  the  tropo¬ 
sphere  than  wave  1. 

In  contrast  to  Au(k)  Fig.  26  shows  that  Av(k)  only  has  one  main 
center  which  Is  centered  in  the  polar  regions.  Both  Au(2)  and  Av(2) 
are  shifted  slightly  south  of  Au(l)  and  Av(l)  respectively.  This 
southward  shift  also  seems  to  hold  true  for  Aw(k)  (Fig.  27),  but  is  not 
as  evident  because  of  the  multiple  cells.  In  addition  to  the  height 
difference  of  wave  2  as  compared  to  wave  1,  we  note  (Fig.  27)  that  wave 
2  has  a  dominant  cell  near  the  tropopause  region  at  about  25°  to  30°N 
(Figs.  27c  and  27d).  Wave  1  also  has  this  maximum  area,  but  It  is  In 
the  equatorial  region  and  does  not  seem  to  be  connected  with  the  strato¬ 
sphere  (Figs.  27a  and  27b). 

From  Fig.  28  we  can  see  the  evolution  of  the  velocity  amplitudes 
from  day  16  through  day  30.  Even  though  Au(l)  Is  the  largest  quantity 
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of  the  two  waves  (Fig.  28a),  It  should  be  remembered  that  Au(l)  is 
mainly  in  the  mesosphere  and  Au(2)  is  dominant  In  the  stratosphere. 
Figure  28b  Indicates  a  slightly  different  picture  In  terms  of  dominant 
wavenumber  In  that  after  day  20  Av(2)  has  the  largest  magnitude.  Aw(k) 
shown  In  Fig.  28c  shows  Cl  starting  with  a  dominant  wave  1  vertical 
velocity  shifting  to  a  dominant  wave  2  on  day  24.  This  Is  viewed  with 
Interest  since  It  occurs  during  the  middle  of  the  warming,  but  Its 
total  significance  Is  not  understood.  The  relatively  large  vertical 
velocity  located  at  about  35°N  and  13  Km  seems  to  Indicate  that  wave 
2  is  the  dominant  force  of  the  model.  This  would  seem  only  logical 
since  our  orographic  effects  also  have  a  wave  2  form  (except  that  In 
C2  we  see  waves  1,  3,  and  4  also  contribute  to  the  vertical  velocity 
even  though  it  also  has  the  same  forcing). 

C2  is  similar  to  Cl  In  that  Au(k)  (Fig.  29)  has  three  maxima, 

Av(k)  (Fig.  30)  has  one  maximum,  and  these  maxima  are  in  approximately 
the  same  latitudinal  location.  They  are  also  similar  in  that  wave  1 
tends  to  be  at  a  higher  altitude  than  wave  2,  but  as  the  dominant  wave 
shifts  from  wave  2  to  wave  1,  wave  2  shifts  downward.  A  difference 
between  the  cases  Is  that  Aw( k )  (Fig.  31)  has  two  maxima  showing.  This 
may  be  due  to  the  stronger  vertical  velocities.  A  second  major  differ¬ 
ence  between  *-he  cases  Is  the  pronounced  shift  of  dominant  wavenumber 
from  wave  2  to  wave  1  In  C2. 

The  magnitude  evolution  graphs  show  a  significant  change  in  Au(k) 
(Fig.  32a)  and  Av(k)  (Fig.  32b)  from  dominant  wave  2  to  dominant  wave 
1.  It  Is  Interesting  to  note  that  this  occurs  during  the  period  of  a 
strong  warming,  but  may  not  be  significant  since  the  second  warming 
pulse  of  C2  does  not  have  this  characteristic.  The  decrease  of  wave  2 
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with  the  Increase  of  wave  1  strongly  suggests  a  non-linear  Interaction 
between  the  waves.  This  also  occurs  In  Ay^k)  (Fig.  33a  and  33b).  As 
waves  1  and  3  decline,  waves  2  and  4  Increase  from  day  22  to  day  24. 
Then  wave  2  decreases  rapidly  as  wave  4  continues  to  increase  and  wave 

3  oscillates  upward.  These  types  of  interactions  continue  through  the 
warming  period,  forming  a  complex  pattern  of  non-linear  Interaction. 
The  C2  Aw(k)  Is  one  of  the  few  values  that  have  significant  wave  3  and 

4  values.  It  Is  also  Interesting  to  note  that  the  dominant  vertical 
wave  at  the  tropopause  Is  wave  3  in  contrast  to  Cl  where  it  was  wave 
2. 


6.2  Cross  Spectral  Analysis 
of  u,  v,  w,  T  and 

The  first  quantity  to  be  examined  is  the  zonal  mean  of  the  meri¬ 
dional  transport  of  zonal  momentum  Cuv(k).  From  Fig.  34a  and  34c  it  is 
easy  to  see  that  the  momentum  flux  Is  consistent  with  the  Cl  velocities 
in  that  wave  1  is  dominant  in  the  mesosphere  and  wave  2  is  dominant  in 
the  stratosphere.  In  addition,  wave  2  has  a  significant  flux  of 
momentum  in  the  equatorial  latitudes.  3y  day  24  of  Cl  there  Is  still 
a  strong  positive  flux  of  Cuv(k)  in  the  mesosphere  at  about  55°N,  but 
an  area  of  divergent  momentum  has  developed  at  75°N  in  the  lower  meso¬ 
sphere.  This  same  phenomena  is  also  true  with  wave  2  only  at  a  lower 
altitude  such  that  there  Is  a  significant  area  of  divergenct  momentum 
In  the  troposphere. 

C2  (Fig.  35)  Is  very  different  from  Cl  concerning  Cuv(k).  First, 
waves  1  and  2  are  closer  to  the  same  altitude.  Second,  wave  1  is  pri¬ 
marily  divergent.  Probably  the  most  distinguishing  difference  between 
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the  Cl  Cuv(k)  and  the  C2  C,JV(k)  Is  that  Its  magnitudes  are  smaller  than 
Cl  before  the  warming,  but  significantly  larger  after  some  warming. 
Additionally  C2  has  no  significant  convergence  or  divergence  of  Cuv(k). 
By  comparing  Figs.  25a,  26a,  and  34a,  It  can  be  seen  that  the  maximum 
Cuv(k)  occurs  between  the  maximum  of  Au(k)  and  Av(k).  This  Implies 
that  the  location  and  amplitude  of  Cuv(k)  Is  a  function  of  the  phase 
relation  between  Au(k)  and  A v( k ) .  This  implication  helps  explain  the 
pulsing  of  the  meridional  circulation.  The  time  evolution  of  Cuv(k) 
for  Cl  (Fig.  36a)  shows  a  decline  of  wave  1  throughout  the  critical 
simulation  period  with  wave  2  declining  from  day  16  to  day  22,  then 
Increasing  rapidly  during  the  peak  of  the  warming.  But,  characteristic 
of  Cl,  both  wave  1  and  wave  2  form  divergent  areas  of  Cuv(k)  during  the 
wind  reversal  period.  This  divergent  momentum  Is  consistant  with  the 
findings  of  O'Neill  and  Taylor  (1979). 

Figure  36b  Is  the  C2  maximum  magnitude  evolution  of  Cuv(k).  As  Is 
characteristic  of  C2,  it  depicts  the  more  rapid  fluxuatlons  of  momentum 
transport.  The  two  peaks  of  divergence  of  wave  1  have  a  direct  correla¬ 
tion  to  the  2  warming  periods  that  C2  experienced.  From  Eq.  (58)  (Kao, 
1980) 

_  m 

3  -  V  1  1  i_  [v(k,t)u(-k,t)cos$]  +  !_  [w(k,t)u(-k,t)] 

3t  ‘  a  cos  $  3$  3z 

.  tan_j_  u(k.t)V(-k,t)  }  +  fv  +  C  (0,t)  (65) 

a  1 

expect  that  as  the  divergence  increases,  the  mean  zonal  velocity  would 
Increase.  In  Eq.  (58)  u(k,t),  v(k,t)  and  w(k,t)  one  Fourier  longitude 
transform  of  u,  v  and  w  respectively.  At  first  glance,  this  seems  to 


be  questionable  for  both  Cl  and  C2.  In  Cl  our  largest  divergent  magni¬ 
tude  occurs  on  day  24,  the  time  of  maximum  zonal  velocity  decline.  In 
C2  we  have  maximum  divergence  on  days  30  and  36.  Day  36  is  during  a 
period  of  strong  zonal  mean  decrease  and  day  30  is  a  time  of  changing 
de-acceleration,  but  still  a  period  of  declining  velocities.  However, 
we  should  note  that  the  first  term  of  the  equation  Is  a  gradient  and 
thus  the  equation  does  hold  true  at  about  80°N. 

The  cross-spectrum  of  the  vertical  transport  of  geopotential 
[Cw$(k)]  Cl  shows  (Fig.  37)  most  flux  to  be  north  of  60°N.  Again  It 
should  be  noted  that  even  though  wave  1  has  over  twice  the  magnitude  as 
wave  2  on  day  18  (Figs.  37a  and  37c),  It  is  wave  2  that  Is  In  the  area 
of  the  wind  reversal.  8y  day  24,  which  Is  the  middle  of  the  major 
warming,  these  distinct  maximum  areas  of  wave  1  and  wave  2  (Figs.  37b 
and  37d)  have  become  much  less  organized,  with  wave  1  propagating  Into 
the  stratosphere. 

In  C2,  the  CW(^{ k )  does  not  have  the  degree  of  altitude  separation 
of  Cl,  but  the  wave  2  magnitude  is  much  larger,  and  thus  could  be  over¬ 
shadowing  the  altitude  differences  of  wave  1  and  wave  2. 

By  day  30  of  C2  (Fig.  38b)  and  Cl  (Fig.  37b)  there  Is  an  area  of 
negative  Cw<j>(  1)  which  has  shifted  downward  slightly  into  the  strato¬ 
sphere.  As  we  look  at  Fig.  39,  these  areas  of  negative  CW(J)(1)  take  on 
an  increased  significance.  This  Is  due  to  their  direct  correlations  to 
the  PNJ  decreases  In  both  Cl  and  C2. 

Figure  39  Indicates  a  major  difference  between  Cl  and  C2  concern¬ 
ing  CW(j,(k).  In  Cl  both  wave  1  and  wave  2  have  a  general  decrease  in 
magnitude  (Fig.  39a).  But  C2  shows  (Fig.  39b)  a  non-linear  interaction 
between  the  two  waves  from  day  24  to  day  32. 


Figure  40  shows  that  the  areas  of  maximum  zonal  mean  meridional 
transport  of  sensible  heat  have  very  similar  locations  to  those  of 
Cvu(k)  (Fig.  34)  and  Cw^,(k)  (Fig.  37).  It  Is  also  Interesting  to 
note  that  for  Cl  all  three  fluxes  have  shown  a  downward  shift  of  the 
flux  with  time,  particularly  CVu(2)  and  Cvj(2). 

In  contrast,  C2  shows  a  slight  upward  shift  of  the  fluxes  (Figs.  35 
and  38)  with  time,  especially  Cv-j-(k)  (Fig.  41). 

The  evolution  of  Cvj(k)  forms  a  different  configuration  than 
C vu ( k )  and  Cw^( k ) .  However,  In  Cl  there  Is  a  good  correlation  between 
the  increase  of  Cvj(2)  (Fig.  42a)  and  the  poleward  shift  of  the  maximum 
temperature  changes  (Fig.  16c).  The  Cvy(k)  of  C2  (Fig.  42b)  seems 
unique  In  that  there  were  no  good  correlations  evident,  with  the 
exception  of  the  non-linear  interaction  of  waves  1  and  2. 

6.3  Polar  Projections  of  Waves 

The  advantage  of  the  polar  projection  of  the  v'  Is  that  we  can  see 
the  horizontal  distribution  and  the  perturbation  velocities  at  different 
levels  in  the  model.  By  counting  the  number  of  alternating  negative 
and  positive  cells  of  meridional  perturbations,  we  can  Indirectly  show 
the  number  of  waves  present. 

We  have  chosen  7.5  Km,  40.5  Km,  and  67.5  Km  to  represent  the  tropo¬ 
sphere,  stratosphere,  and  mesosphere  respectively.  For  the  troposphere, 
we  selected  three  days  to  show  the  wave  evolution,  but  since  the  upper 
atmosphere  has  less  wave  fluxation  we  selected  only  two  days. 

At  7.5  Km  In  Cl,  day  14  (Fig.  43a)  we  have  wave  2  nearly  symmet¬ 
rical  positioned  around  the  North  Pole.  By  day  18  (Fig.  43b),  the  wave 
amplitudes  have  nearly  doubled  with  a  weak  wave  three  forming  near  0® 


longitude.  By  day  30  (Fig.  43c),  the  third  wave  has  strengthened  and 
the  North  Pole  Is  under  the  Influence  of  a  weak  ridge.  In  the  strato¬ 
sphere  on  day  13  (Fig.  44a),  there  Is  a  clear  wave  2  with  the  turbulent 
velocities  nearly  equal  to  the  zoan  mean  velocities.  By  day  30  (Fig. 
44c),  the  waves  have  weakened  and  appear  to  possibly  be  In  a  transition 
from  wave  2  to  wave  1.  The  mesosphere  starts  with  a  wave  l  (Fig.  44b) 
and  shifts  to  a  weak  wave  2  by  day  30  (Fig.  44d).  The  dominant  wave 
has  shifted  Into  the  Eastern  Hemisphere. 

C2  shows  a  similar  wave  configuration  on  day  14  at  7.5  Km  (Fig. 
45a).  This  pattern  Is  still  similar  on  day  18  though  less  symmetrical 
than  Cl  (Fig.  45b).  By  day  30,  however,  the  troposphere  has  shifted  to 
a  wave  4  configuration  (Fig.  45c)  Instead  of  the  wave  3  pattern  of  Cl. 
At  40.5  Km  on  day  18  (Fig.  46a),  there  Is  a  symmetrical  wave  2. 

However,  by  day  30  (Fig.  46c),  these  waves  have  begun  to  shift  to 
a  wave  1  centered  at  the  North  Pole.  The  mesosphere  starts  with 
essentially  a  wave  1  configuration  on  day  18  (Fig.  46b).  The  wave  1 
continues  to  dominate  the  67.5  Km  level  through  the  warming  period  even 
though  It  Is  distorted  some  by  a  weak  wave  2  (Fig.  46d). 


CHAPTER  7 


CONCLUSION 

Using  an  orographic  forcing  of  a  combination  of  wavenumbers  1  and 
2  In  the  troposphere,  major  (case  1)  and  minor  (case  2)  stratospheric 
warmings  were  simulated  with  a  31-layer  primitive  equation,  spectral 
model  consisting  of  a  5-layer  troposphere,  12-layer  stratosphere  and 
14-layer  mesosphere.  An  analysis  of  the  data  resulting  from  the  simu¬ 
lations  Indicated  several  similarities  and  differences  between  the 
warmings.  Case  1  had  a  weak  polar  vortex  and  case  2  had  a  strong 
polar  vortex.  The  effect  of  the  forcing  on  the  zonal  winds  was  to 
Induce  wave  motion,  meridional  and  vertical  flux  of  momentum  and  sensl 
ble  heat  In  the  model  atmosphere.  The  vertical  flux  of  geopotentlal 
began  near  45°N,  then  shifted  northward.  On  about  an  eight  to  twelve 
day  cycle,  the  meridional  and  vertical  fluxes  of  momentum  and  sensible 
heat  reach  a  peak  In  magnitude  that  coincided  with  the  peak  of  the 
meridional  circulation.  With  the  first  peak  of  the  meridional  clrcula 
tlon,  the  polar  night  jet  became  a  distinct  feature.  The  second  peak 
of  meridional  circulation  coincided  with  the  wind  reversal  In  case  1 
and  a  significant  decrease  In  the  mean  polar  winds  In  case  2.  The 
meridional  circulation  did  not  reach  a  third  peak  In  case  l,  but  In 
case  2,  It  was  associated  with  a  further  de-acceleration  of  the  polar 
night  jet. 

The  similarity  of  the  configurations  of  the  amplitude  maxima 
evolution  of  the  turbulent  transfer  terms  Is  evident  In  both  cases. 
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The  reason  for  this  Is  probably  a  result  of  the  Indirect  relationship 
between  the  meridional  flux  of  sensible  heat  and  the  meridional  trans¬ 
port  of  zonal  momentum  that  Is  shown  to  exist  as  a  result  of  the 
theorems  of  non-interaction  and  the  continuity  equation.  These  rela¬ 
tionships  would  also  provide  a  possible  explanation  for  the  periodicity 
of  the  Indirect  merlodlonal  cell  circulation,  as  well  as  the  semi- 
11  near  temperature  changes  In  the  stratosphere. 

In  both  cases  we  found  that  wave  number  1  waves  were  dominant  In 
the  mesosphere  and  wave  number  2  waves  were  dominant  In  the  strato¬ 
sphere.  However,  during  the  warming,  wave  number  1  waves  became 
significant  In  the  stratosphere.  Case  2  exhibited  more  evidence  of 
nonlinear  Interaction  than  case  1. 

Ts*  -alldlty  of  the  assumptions  made  In  deriving  the  relationship 
betv»fc  *.,e  vertical  flux  of  geopotential  and  the  meridional  flux  of 
sensible  heat  from  the  linearized  governing  equations  was  tested  with 
the  use  of  the  data  generated  from  the  numerical  model  analyzed.  We 
found  that  for  the  quasi -stationary  state  there  was  some  correlation 
between  the  linearized  theory  and  the  zonal  mean  latitude-height 
positions  of  the  fluxes.  However,  for  the  nonstationary  condition  of 
the  warming  process,  the  correlation  become  poor.  This  tends  to  raise 
questions  about  the  validity  of  the  application  of  the  relationship 
between  the  vertical  flux  of  geopotential  and  meridional  flux  of 
sensible  heat  derived  from  the  llnarlzed  governing  equations  to  non¬ 
linear  processes,  and  therefore,  the  concept  of  the  critical  level  to 
the  process  in  the  developing  stage  of  strong  stratospheric  warmings. 
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Fig.  2.  Latitude-height  sections  of  Initial  mean  zonal  wind 

(x  10  m  s"M  for  (a)  case  1;  and  (b)  case  2  (Koermer,  1980). 

Used  by  permission. 
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Fig.  3.  Maximum  longitudinal  surface  height  distribution  at  45*N 
(Koermer,  1980).  Used  by  permission. 


Latitude-height  sections  of  the  zonal  mean  variance  of  u 
for  (a)  case  1,  day  18;  (b)  case  1,  day  30;  (c)  case  2, 


-uJ. 


23  3J  46  50  1  60  ‘  70“  >0  90 

LATITUDE  iOEG) 


1  1  ■  ~  ~i~.  ■  ,  ■  ,  ■  ■  :  ■  1 

®0  10  20  10  40  so  60  70  80  90 

LATITUDE  IDECI  (d) 


Latitude-height  sections  of  the  zonal  mean  zonal  velocity 
(x  10  m  s"1)  for  (a)  case  1,  day  18;  (b)  case  1,  day  30; 
(c)  case  2,  day  26;  (d)  case  2,  day  30. 
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g.  7.  Latitude-height  sections  of  the  zonal  mean  variance  of  the 
meridional  wind  for  (a)  case  1,  day  18;  (b)  case  1,  day 
30;  (c)  case  2,  day  26;  (d)  case  2,  day  30.  Units  for 
(a)  and  (b)  are  (x  10?  m2  s"2).  Units  for  (c)  and  (d)  are 
(x  103  m2  s-2). 
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Fig.  10.  (a)  and  (b)  are  latitude-height  sections  for  case  1  with 

dashed/sol  id  lines  being  zonal  mean  meridional  velocity 
(m  s"M  and  zonal  mean  vertical  velocity  (x  10~2  m  s"1) 
respectively,  (c)  case  1  zonal  mean  meridional  velocity 
maximum  evolution;  (d)  (Koermer,  1980)  case  1  latitude¬ 
time  section  of  the  zonal  mean  meridional  velocity  at 
40.5  Km.  Units  for  (c)  and  (d)  are  (m  s-1).  Used  by 
permission. 
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Case  1  mean  vertical  velocities  ( 10~2  m  s-l)  (a)  velocity 
maxima  evolution;  (b)  (Koermer,  1980)  latitude-time  section 
at  40.5  km.  Used  by  permission. 


Fig.  12.  Same  as  Fig.  10  except  for  case  2 


Imposed  on  zonal  mean  meridional  flux  of  zonal  momentum  (solid  line).  Units  are  (m  s-1) 
and  (m^  s'*)  respectively;  (c)  zonal  mean  meridional  velocity  (dashed  line)  superimposed 
on  the  zonal  mean  vertical  transport  of  zonal  momentum  (solid  line).  Units  are  (m  s'1) 


Fig.  15.  Evolution  of  mean  zonal  wind  maxima  (m  s"*)  with  solid  line 
being  polar  night  jet,  dash  line  being  mesospheric  jet, 
dotted  line  being  tropspherlc  jet,  and  dot  dashed  line 
being  mid-latitude  easterlies  for  (a)  case  1;  (b)  case  2. 
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Fig.  16.  Temperature  evolution  charts  (solid  line  Is  stratosphere 
warming,  dashed  line  Is  mesosphere  cooling)  for  case  1  of 
(a)  temperature  maximum  change  with  time  (°C);  (b) 
height-time  evolution  of  temperature  maximum  changes;  (c) 
latitude-time  evolution  of  temperature  maximum  changes. 


Fig.  18.  Time  charts  of  the  zonal  mean  vertical  flux  of  geopoten¬ 
tial  maxima  for  case  1  depicting  (a)  evolution  of  maxima 
{x  10z  m3  s"3);  and  (b)  latitude  location  of  the  maxima. 
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Fig.  20. 


Time  charts  of  the  zonal  mean  meridional  flux  of  zonal 
momentum  for  case  1  depicting  (a)  evolution  of  maxima 
(x  102  m2  s"2h  and  (b)  latitude  location  of  the  maxima. 


Fig.  22.  Magnitude-time  charts  of  the  zonal  mean  of  the  vertical 
flux  of  zonal  momentum  maxima  (m2  s"2)  for  (a)  case  1 
where  solid  line  Is  in  the  upper  stratosphere  from  60°N  - 
80°N,  dashed  line  Is  in  low  mesosphere  from  37°N  -  47°N; 
and  (b)  case  2  where  solid  line  is  in  the  upper  stratopause 
lower  mesopause  from  65°N  -  838N,  dashed  line  is  In  the 
lower  mesosphere  from  44°N  -  65°N  dotted  line  is  at  48°N 
and  70  Km. 
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27.  Latitude-height  sections  of  A^k)  (x  10"2  m  s"1)  for  case 
1  (a)  day  18,  wave  1;  (b)  day  24,  wave  1;  (c)  day  18,  wave 
2;  and  (d)  day  24,  wave  2. 


Fig.  29.  Latitude-height  sections  of  A„(k)  for  case  2  (a)  day  26, 
wave  1;  (b)  day  30,  wave  1;  (c)  day  26,  wave  2;  and 
(d)  day  30,  wave  2.  Units  are  m  s"l. 
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Fig.  30.  Latitude-height  sections  of  AVU)  for  case  2  (a)  day  26, 
wave  1;  (b)  day  30,  wave  1;  (c)  day  26,  wave  2;  (d)  day 
30,  wave  2.  Units  are  m  s_1. 
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Latitude-height  sections  of  /^(Ic)  (x  10"2  m  s-1)  for  case 
2  (a)  day  26,  wave  1;  (b)  day  30,  wave  1;  (c)  day  26, 
wave  2;  and  (d)  day  30,  wave  2. 


Latitude-height  sections  of  Cyv ( k )  (x  10  m £  s"^)  for  case  1 
(a)  day  18,  wave  1;  (b)  day  24,  wave  1;  (c)  day  18,  wave  2; 
and  (d)  day  24  wave  2. 


Latitude-height  sections  of  Cyv ( k )  (x  10  s"^)  for  case  2 

(a)  day  26,  wave  1;  (b)  day  3u,  wave  1;  (c)  day  26,  wave 
2;  and  (d)  dfiy  30,  wave  2. 
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Fig.  44.  Polar  stereographic  projections  of  meridional  velocity 

deviations  (x  10  m  s’M  for  case  1  on  (a)  day  18,  40.5  km; 
(b)  day  18,  €7.5  km;  (c)  day  30,  40.5  km;  and  (d)  day  30, 
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Polar  *t«r*oqrapMc  projections  of  aerlodlonal  velocity 
deviations  U  10  m  *_1)  for  case  2  on  (a)  day  18,  40.5  kn 
(t>)  day  18,  67.5  k«;  (c)  day  30,  40.5  ka;  and  (d)  day  30, 
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